Renormalization of the elementary excitations in hole- and electron-doped cuprates 

due to spin fluctuations 
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Extending our previous studies we present results for the doping-, momentum-, frequency-, and 
temperature- dependence of the kink-like change of the quasiparticle velocity resulting from the 
coupling to spin fluctuations. In the nodal direction a kink is found in both the normal and super- 
conducting state while in the antinodal direction a kink occurs only below T c due to the opening of 
the superconducting gap. A pronounced kink is obtained only for hole-doped, but not for electron- 
doped cuprates and is characteristically different from what is expected due to electron-phonon in- 
teraction. We further demonstrate that the kink structure is intimately connected to the resonance 
peak seen in inelastic neutron scattering. Our results suggest similar effects in other unconventional 
superconductors like S^RuCU. 
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I. INTRODUCTION 

Elementary excitations in the cuprates are of central 
interest in order to learn more about the correlations and 
the pairing mechanism for superconductivity. For exam- 
ple, it is well-known that the understanding of the ele- 
mentary excitations in conventional superconductors like 
lead as measured by tunneling spectroscopy played the 
crucial role in accepting the picture of electron-phonon- 
mediated Cooper-pairing 1-3 . In the high-T c cuprates one 
expects that due to the presence of antiferromagnetic 
spin fluctuations a strong renormalization of the spectral 
density and the corresponding energy dispersion may oc- 
cur. It was shown by several groups 4,5 that the so-called 
'dip-hump' structure seen in tunneling experiments 6 can 
be explained in terms of the feedback effect of supercon- 
ductivity arising from the structure in the gap function 
A(k, ui). Moreover, it has been argued that this struc- 
ture reflects the effective pairing interaction and points 
towards a spin-fluctuation-mediated Cooper-pairing. Re- 
cent developments in angle-resolved photoemission spec- 
troscopy (ARPES) allow to study the elementary exci- 
tations directly and in much more detail. In particular, 
the anisotropy of the elementary excitations close to the 
Fermi energy in the different parts of the Brillouin Zone 
(BZ) are studied. This is important, since the coupling of 
the quasiparticles to spin fluctuations varies at different 
parts of the BZ. Furthermore, the understanding of the 
structures seen by ARPES and their doping, momentum, 
and temperature dependence will help to understand the 
role played by spin fluctuations in contrast to phonons re- 
garding the formation of superconductivity in the high-T c 
cuprates. 

1. Theoretical background: Similar to the electron- 
phonon case, the coupling between the quasiparticles and 
the spin excitations should influence characteristically 
the energies Wk of the hole- or electron carriers 



Here, €k refers to the bare energy dispersion of the quasi- 
particles assuming no interaction with the spin excita- 
tions or phonons, i.e. a tight-binding energy disper- 
sion for the Cu02-plane. In general, the self-energy 
S(k, oj) results from the coupling of the particles to spin 
excitations, see Fig. 1 for an illustration. Obviously, 
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FIG. 1. Illustration of the coupling between holes or elec- 
trons and spin fluctuations characterized by the susceptibility 
x(q, w)- This leads to Dyson's equation G -1 = Gq 1 — E de- 
scribing the relation between the bare Green's function Go 
and the renormalized one G(k, ui). U (f7 e ff) denotes an effec- 
tive (renormalized) coupling constant. 

the self-energy is a functional of the spin susceptibil- 
ity, £ = 111 its simplest form, the latter is of 
Ornstein-Zernicke form 7 that allows for a sharp enhance- 
ment of fluctuations near the antiferromagnetic wave vec- 
tor Q= (it, 7r): 
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Here, xq is the value of the static spin susceptibility at 
the wave vector Q, £ is the magnetic correlation length, 
and u s f is the characteristic frequency of spin fluctua- 
tions (roughly the peak position in the imaginary part of 
Eq. (2)). Due to the fact that x(q, u>) has only structures 
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FIG. 2. Illustration of the anisotropy of the elementary ex- 
citations. The solid line denotes the calculated Fermi surface 
for the hole-doped cuprates in the first BZ using Eq. (3). 
The dashed lines refer to the magnetic Brillouin Zone that 
crosses the electronic Fermi surface close to the (tt, 0) points 
exactly at the 'hot spots' where the antiferromagnetic wave 
vector Q = connects two pieces of the Fermi surface. 

At the (tt, 0) point and along the diagonal the wave vector 
Q connects quasiparticle states below the Fermi level only. 
This allows us to define three characteristic regions 1,2,3 at 
the Fermi level. 

around u> s f and wave vector Q, their influence on the el- 
ementary excitations is expected to be very anisotropic 
at different parts of the Brillouin Zone. 

In order to illustrate the anisotropy of the elementary 
excitations we show in Fig. 2 the calculated Fermi surface 
for hole-doped cuprates 8 . For the calculation we use the 
tight-binding energy dispersion for the CuCVplanc 

ek = — 2t (cos k x + cos k y ) + At' cos k x cosfc y — \i . (3) 

Here, t and t' refer to the hopping of a hole (electron) 
between nearest, next-nearest sites on the square lattice, 
and /j, is the chemical potential that defines the doping. 
In Eq. (3) and in the following we set the lattice constant 
to unity. One clearly sees that the scattering of quasi- 
particles by spin fluctuations in cuprates is anisotropic. 
First, the antiferromagnetic wave vector Q = (w, tt) con- 
nects exactly quasiparticles at the Fermi level close to the 
(tt, 0) points of the first Brillouin Zone. These quasipar- 
ticles experience the strongest coupling to antiferromag- 
netic spin fluctuations. Quasiparticles at the diagonals 
of the BZ are not connected by the Q and thus have 
smaller scattering by spin fluctuations at the Fermi level. 
The corresponding point on the Fermi surface is called 
cold spot or nodal point. Furthermore, the quasiparticle 
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FIG. 3. Calculated bare tight-binding energy dispersion for 
the hole-doped cuprates in the normal state using Eq. (3). 
The dashed curve illustrates the changes due to renormaliza- 
tion, (j k = e k + ReS(k, u>). Due to structure in ReE the en- 
ergy dispersion shows a kink feature along the (0, 0) — > (tt, n) 
direction in the first Brillouin Zone. 

states at the Fermi level close to the M point in the first 
BZ are usually called antinodal points. Note that the 
quasiparticle states at cold and hot spots which are con- 
nected by Q lie close to the Fermi level. This will be 
important later for discussing the kink feature. There- 
fore, we may define three characteristic regions at the 
Fermi surface regarding their sensitivity to coupling to 
antiferromagnetic spin fluctuations. 

2. Experimental findings: Let us now discuss the 
experimental situation close to the Fermi level as mea- 
sured by ARPES. The combined study of energy distri- 
bution curves (EDC) and momentum distribution curves 
(MDC) allows the study of the quasiparticle excitations 
close and well below the Fermi level up the energies of 
200 meV. Most importantly, recent ARPES studies by 
various groups 10-13 reveal a kink structure in the energy 
dispersion at energies about 60 ± 15meV below the Fermi 
energy and along the nodal direction (0, 0) — ► (7T, 7t) di- 
rection of the first BZ. Furthermore, this kink structure 
appears in the normal state and almost does not change 
if one enters the superconducting state. 

We illustrate the formation of the kink in Fig. 3 where 
we show the unrenormalizcd tight-binding energy dis- 
persion £k in the normal state along the route (0,0) — > 
(tt,0) (tt.tt) (0,0) of the first BZ. Below E f , the 
measured energy dispersion Wk along (0, 0) — > (tt, tt) has 
changed compared to the bare tight-binding case due to 
the self-energy corrections as it follows from Eq. (1). 
We define the kink at the inflection point of the dashed 
curve where the renormalized dispersion tends to ap- 
proach again the bare dispersion ek- Recent studies by 
Dessau and co-workers 14 reveal no kink structure along 
(tt, 0) — ► (tt, tt) direction in the normal state, but only in 
the superconducting state. 
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3. Doping and k-dependence of the kink: Pre- 
viously we have shown 15 that from the momentum and 
energy conservation one expects changes of the quasi- 
particle velocities = de/dk in the nodal direction 
for u ~ oj s f and k k' = q ~ Q. This leads to 
the kink in the (0, 0) — > (n, tt) direction that has been 
observed by angle-resolved photoemission spectroscopy 
(ARPES) 10-13 . Here, we extend our previous analy- 
sis and discuss in more detail the doping dependence 
and the anisotropy of the kink feature as a function of 
(k — k F ) and energy w k in the normal and supercon- 
ducting state, for different routes k kp in the Brillouin 
zone: (0,0) -> (tt.tt), (0,0) (0,tt), and (tt,0) -> 
Due to the fact that the superconducting order param- 
eter has nodes along (0,0) — > (ir,n) and maxima for 
(0, 0) — > (0, 7r) one expects additional changes in the 
kink structure as observed recently 16 ' 14 . This should help 
to contrast renormalization due to spin fluctuations to 
what we get due to electron-phonon interaction. Fur- 
thermore, we investigate the interdependence of the ele- 
mentary excitations with the resonance peak that is ob- 
served in hole-doped cuprates by inelastic neutron scat- 
tering (INS) below T c . It follows also from Eq. (1) and 
Im£(k, uS) oc lu 2 — > u> that the temperature dependence 
of the kink feature should reflect characteristically the 
coupling to spin fluctuations and be different from the 
electron-phonon coupling case. 

Note that in the case of electron-doped cuprates nei- 
ther a resonance peak nor a kink feature is present. We 
will also show that the kink feature is not restricted to 
the cuprates, but is also expected for other novel super- 
conductors where quasiparticles couple strongly to spin 
excitations. A particularly interesting case might be 
Sr2RuC>4 with large anisotropic behavior of X zz (q, w) and 
in-plane xh (q, w) 17 ' 18 . Generally, a strong nesting be- 
havior of the Fermi surface might yield pronounced kink 
features. 

This paper is organized as follow: In Section II we 
present the theory, and in Section III we discuss our re- 
sults for the kink structure, its doping- and tempera- 
ture dependence for hole- and electron-doped cuprates. 
Finally, in Section IV we summarize our analysis and 
contrast renormalization due to spin fluctuations and 
phonons. 

II. THEORY 

A. One-band Hubbard model: The elementary ex- 
citations and the spin excitations are key quantities de- 
termining the properties of the cuprates and other su- 
perconductors with strong correlations and magnetic ac- 
tivity. The quasiparticles, holes or electrons, interact 
strongly with spin fluctuations and also with phonons. 
However, phonons and spin fluctuations differ with re- 
spect to their doping dependence and anisotropy. This is 
clearly demonstrated by neutron scattering experiments, 



for example. Also the different behavior of hole- and 
electron-doped cuprates and the feedback of supercon- 
ductivity on x(q, uj) is important. 

In this paper we employ an effective one-band Hamil- 
tonian. This is justified because upon hole doping an- 
tiferromagnetism disappears due to Zhang-Rice singlet 
formation and quenching of Cu-spins. In this one-band 
picture the Coulomb interaction between the quasiparti- 
cles refers to an effective interaction (i.e. the Hubbard 
U) within the conduction band. Then, further doping in- 
creases the carrier mobility and a system of strongly cor- 
related quasiparticles occurs. In the overdoped case less 
magnetic activity is present yielding usual Fermi liquid. 
We assume U ~ W/2 (W — bandwidth) independent of 
the doping concentration. 

The main physics of a single Cu02-plane is the two- 
dimensional one-band Hubbard model given by 

H = - tl J ( c tr c 3° + c % c ™) +U Y1 ni T n a - ( 4 ) 

(ij) a i 

where c+ creates an electron with spin a on site i, U 
denotes the on-site Coulomb interaction, and Uj is the 
hopping integral. After diagonalization of the first term, 
one arrives at the bare tight-binding energy dispersion 
given by Eq. (3). The description of the electron- and 
hole-doped cuprates within a one-band approximation 
is possible if one takes into account different parame- 
ters and quasiparticle dispersion 19 . Note that in the 
case of electron doping the electrons occupy the cop- 
per ci-band, while in the hole-doped case holes refer 
mainly to the oxygen p-states yielding different disper- 
sion parameters. Furthermore, the energy dispersions 
for optimally hole-doped La2- :r Sr :r Cu04 (LSCO) and 
electron-doped Nd2- £! ;Ce a ;Cu04 (NCCO) behave differ- 
ently around (tt, 0) point. While in the case of LSCO 
the flat band (leading to the van Hove singularity in 
the density of states) lies close to the Fermi level, in 
NCCO the flat band is approximately 300 meV below 
the Fermi level. Then, using t = 250meV and t' = O.li, 
one describes the hole-doped LSCO dispersion, whereas 
t = 138meV and t' = 0.3i are needed for the descrip- 
tion of the electron-doped NCCO compound fitting ear- 
lier photoemission data. These parameters including an 
intermediate coupling U = 4i, will be used for our cal- 
culations of various physical quantities in the normal 
and superconducting state of electron- and hole-doped 
cuprates 20 . 

B. Generalized Eliashberg equations: In this one- 
band model, we assume that the same electrons (holes) 
are participating in the formation of antiferromagnetic 
fluctuations and in Cooper-pairing due to the exchange 
of these fluctuations. Thus, both the magnetic suscep- 
tibility and the quasiparticle self-energy must be cal- 
culated self-consistently. This is possible in the FLEX 
approximation 21-24 . In this approach the dressed one- 
electron Green's function are used to calculate the charge 
and spin susceptibilities. These susceptibilities are then 
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used to construct an effective Berk-Schrieffer-like 25 pair- 
ing interaction V e s describing the exchange of charge and 
spin fluctuations. The generalized Eliashberg equations 
are derived in Appendix A. In order to demonstrate the 
significant role of V c g in our work, we also show in Fig. 13 
its corresponding diagramatic representation. Note, in 
general, if the Cooper-pairing and the effective pairing 
potential V s are generated by the same quasiparticles 
(solid lines in Fig. 13), strong self-energy and feedback 
effects on G(k, uj) and %(q, Q) are expected 26 . 

To be more precise, we write down the quasiparticle 
self-energy components Y, v (v = 0, 3, 1) with respect to 
the Pauli matrices t„ in the Nambu representation 28 ' 29 , 
i.e. So = w(l — Z) (mass renormalization), S3 = £ 
(energy shift), and Si = <j> (gap parameter). They are 
given by 

/•OO 

£„(k, w) = N- 1 Y. dil V cS (k - k', fi) 
k' Jo 

/ + OC 
dw'J(w,n,w')A I ,(k / ,w / ) (5) 
-00 

with 

V cS = [P s (k - k', SI) - (5 vl - 5 v0 - 6 v3 ) P c (k - k', 0)] . 

(6) 

This is a generalization of figure 1. P s and P c denote 
the spectral density of the spin and charge excitations, 
respectively, and arc defined in Eqs. (A3) and (A4). 
The second part of Eq. (5) is given in Eqs. (A5) and 
(A6). It is interesting to remark that the above formu- 
lae remain valid even if the elementary excitations and 
the magnetic activity that controls V e e would result from 
different quasiparticles. 

The generalized Eliashberg equations allow us to calcu- 
late all the properties of the system self-consistently like 
superconducting phase diagram, elementary excitations, 
the superconducting order parameter, energy dispersion 
and dynamical spin susceptibility, for example 30 ' 9 . Note, 
in Eq. (5) self-energy effects due to phonons are ne- 
glected. Their contribution will be discussed in section 
IV and in Appendix B. 



III. RESULTS AND DISCUSSION 

A. Anisotropic Renormalization 

1. Nodal direction: We start the discussion analyz- 
ing the spectral density of hole-doped superconductors 
in the normal state. The spectral density reveals the el- 
ementary excitations and in particular the renormalized 
energy dispersion. First, we present our results for the 
spectral density along the nodal (0, 0) — ► (n, ir) direction 
in the first BZ. 




Energy (meV) 



FIG. 4. Calculated spectral density N(k,Lj) in the normal 
state along the nodal (0, 0) — ► (tt, it) direction (from left to 
right) as a function of frequency in the first Brillouin zone 
(BZ). The peak positions (connected by the solid line to guide 
the eye) refer to the renormalized energy dispersion tJk- One 
clearly sees the kink structure at an energy approximately 
•jJki-nk — 75 ± 15meV that results from coupling of the quasi- 
particles to spin fluctuations. 

In Fig. 4 we show the calculated spectral density 
iV(k,u;), i.e. the local density of states, as a function 
of frequency and momentum k — kp. The peak posi- 
tions correspond to the renormalized energy dispersion. 
Due to coupling of holes to antiferromagnetic spin fluc- 
tuations the quasiparticle dispersion changes its slope 
and shows a pronounced kink feature at the energy 
Ukink ~ 75 ± 15meV. 

How can one understand the kink feature in a simple 
way? At the first glance the occurrence of a kink in the 
nodal direction seems to be surprising, since the main 
interaction of the carriers with spin fluctuations occurs 
at the hot spots while the kink feature is present along 
the diagonal of the BZ close to the cold spots. This 
argument, however, considers only the quasiparticles ex- 
actly at the Fermi level. Away from the Fermi level but 
close to it (along (0,0) — ► (tt,7t)) the quasiparticles cou- 
ple strongly to spin fluctuations, as can be seen from 
Fig. 2. Most importantly, as follows from Fig. 2, the 
largest scattering will occur at values of k — kp = Q and 
w = u) s f. To be more precise, let us rewrite Eq. (5) as 

S(k,^„) = -T 2 ^^oG(k-k',z Wll -z^ m )f (7 2 

Wm,^m k',q 

x ^Tr [f G(k + q, iw m + w m )f G(q, iu) m )] (7) 

and approximate the Green's function by its non- 
interacting part 

G(k, tw„) m G (k, tuj m ) = — — — 2 . (8) 

with = e 2 . + 4>1.. Thus, after little algebra one obtains 
on the real axis 32 
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FIG. 5. Spectral density N(k, to) as a function of frequency 
along the (0, 0) — > (0, tt) direction of the first BZ in the normal 
state calculated from the generalized Eliashberg equations. 
Again, the peak positions reveal the renormalized energy dis- 
persion cdk- A kink occurs at similar energy as in the nodal 
direction. Because of inelastic scattering of holes on spin fluc- 
tuations close to (0, 7r), N(k,u) becomes also broader. Note 
that, in contrast to the nodal direction, one does not cross 
the Fermi level in the (0, 0) — > (0, tt) direction. Instead, one 
reaches the flat part of the tight-binding band. 



£(k,w) 
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duj 



,Im x RP A(k-k'X) 
w — u' — £k' 



. Ul \ I u> — U) 

coth — — — tanh — — 

2T J \ 2T 



(9) 



The imaginary part of the spin susceptibility may be 
obtained within the RPA and is approximately given 
by Ornstein-Zernicke expression (see Eq.(2)). This self- 
energy £ now enters Eq. (1). It is important that the 
self-energy is mainly frequency-dependent, while the bare 
dispersion of the carriers is not. Already in the normal 
state, S(k, w) has a maximum reflecting a corresponding 
maximum of Im \ at q w Q and u' ~ u s f. Note that 
u 8 f can be determined according to Moriya and Ueda 
and paramagnon theory from xrpa — Xo/(l — U\o) or 
equivalentlt from x t 7 1 (q — Q,uj ~ cj s f) —U = 33 . Ob- 
viously, uj s f is strongly doping dependent. This will be 
discussed later. Then, the kink position follows from the 
pole of the denominator of Eq. (9). This leads to the 
'kink condition' 



E, 



k-Q 



sf(x) 



(10) 



This gives an estimate of the position of the kink. Fur- 
thermore, since the superconducting gap is zero for w = 
0, but not for to — u> s f, the kink feature along the nodal 
direction (0, 0) — -> (tt, tt) will change only slightly below 
T c . This we have demonstrated previously 15 . 

2. (0, 0) — ► (0, tt) direction: In order to see whether 
the kink feature is present in other directions of the Bril- 
louin Zone, we show in Fig. 5 the evolution of the spectral 
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FIG. 6. Calculated spectral density N(k, lo) along the 
antinodal (tt, 0) —> (tt, tt) direction as a function of frequency 
in the first BZ in the normal (a) and superconducting (b) 
state. Due to the flat band close to the Fermi level the spec- 
tral density shows no kink structure in the normal state. Be- 
low T c the superconducting gap 4>(to) opens yielding a kink 
structure in the spectral density that occurs at the energies 
tOkink ~ 50 ± 10 meV at optimal doping. 



density along the (0, 0) — ► (0, tt) direction. Despite of the 
fact that along this direction we do not cross the Fermi 
level, the kink feature is still present and is found at an 
energy similar to the one for the nodal (0, 0) — ► (tt, tt) 
direction. This indicates that the occurrence of the kink 
feature is not connected to some specific conditions which 
might be present only along the (0, 0) — > (tt, tt) direction. 
Instead, the kink is characteristic for all direction where 
k kp ~ Q and to ~ lo s j. Also below T c we find that 
the kink feature is present in the (0, 0) — > (0, tt) direction 
(not shown). Note, that our results are in fair agreement 
with experimental data 12 . 

3. Antinodal direction: In Fig. 6(a) we show our 
results for N(k,uj) along the (tt, 0) — > (tt,tt) route, i.e. 
the antinodal direction, of the first BZ in the normal 
state. Note, that the spectral density at the (0, tt) point 
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is broader than at the antinodal point due to stronger 
coupling to spin excitations peaked at q = Q = (n, n) 
as discussed in Fig. 2. Clearly, no kink is present. The 
absence of a kink structure can be explained with the 
flat structure of the CuC-2-plane around the M point (see 
Fig. 3). Simply speaking, for a flat band the frequency 
dependence of £ in Eq. (1) does not play a significant 
role and therefore no change of the velocity and no kink 
structure is present. 

What does happen in the superconducting state? Be- 
low T c the superconducting gap <j)(k,ui) opens rapidly 
for decreasing temperature T and becomes maximal in 
momentum space around the M point reflecting the mo- 
mentum dependence of the effective pairing interaction 
(see Eq. (5)). In addition, due to the frequency depen- 
dence of the gap the flat band around M disappears. 

In Fig. 6(b) we show results for iV(k, u>) at a tem- 
perature T = 0.5T C where the superconducting gap has 
opened. A kink structure around uJkink ~ 50 ± 10 meV 
is present reflecting the magnitude of <fi. Hence, in the 
(it, 0) — ► (ir, 7r) direction this kink feature is only present 
below T c and connected to the feedback effect of <p on 
the elementary excitations. We will show later that this 
feedback is also important for the resonance peak seen in 
INS. 

Note, that the superconducting gap </>(k, u>) is calcu- 
lated self-consistently in our theory and reflecting the 
underlying spin fluctuations which dominate the pair- 
ing potential V e s. Therefore, the occurrence of a kink 
structure only below T c in the antinodal direction is a 
direct fingerprint of the spin excitation spectrum 35 . Fur- 
thermore, as we will discuss below, Im x(Q> w ) entering 
in Eq. (9) is peaked at the resonance frequency ui res 
(roughly at u) s f + A). Therefore, the kink condition is 
given by 

U>kink ~ -Efc-Q + ui res (x) . (11) 

In Fig. 7(a) the frequency dependence of Re S(k a , ui) 
in the normal and superconducting state at the antinodal 
point k = k a is shown. Due to the occurrence of the res- 
onance feature in Im x(Q, ui) and the related feedback of 
the superconducting gap <f>(u>), Re S shows a pronounced 
structure below T c at energies of about w res + A . Also 
the corresponding imaginary part, Im S(k = k a ,u)), 
shows a peak below T c (see Fig. 7(b)). This pro- 
nounced behavior is responsible for the kink formation 
along (it, 0) — > (w, n) direction in the BZ. Therefore, 
while the kink features are present along (0, 0) — > (n, n) 
and (-7T, 0) — > (it, tt) directions in the superconducting 
state of hole-doped cuprates, their nature is qualitatively 
different. Along the nodal direction the superconduct- 
ing gap is zero (for ui = 0) and thus the feedback effect 
of superconductivity on the elementary and spin exci- 
tations is small. Therefore, u> s f determines mainly the 
formation of the kink feature. On the other hand, along 
the antinodal direction the gap is maximal and yields a 
strong feedback of superconductivity on \- Thus, in the 



150 




-150 |- - 

-300 -200 -100 100 200 300 

co (meV) 



140 

r-v 120 
> 

2 100 

B 

w 80 









T < Tc 


\. 1 

x = 0.15 


/ 

\ . 


(b) 

antinodal 


\ / T > T c 

\ / 

\. y 



-100 -50 50 100 



co (meV) 

FIG. 7. (a) Calculated 

frequency dependence of Re E(k a ,ai) at the antinodal point 
k a of the first BZ in the normal (solid curve) and supercon- 
ducting state (dashed curve). Due to the feedback effect of 
the superconducting gap 4>(uj), a peak (dip) occurs for ui > 
(uj < 0) which roughly defines the position of the kink struc- 
ture, (b) The corresponding imaginary part at the antinodal 
point Im E(k = k a ,cj) is shown. Again, due to the feedback 
effect of 4>(u>), a maximum occurs below T c . Note, that both 
Re E and Im E are not fully antisymmetric (symmetric) with 
respect to u> at optimum doping x = 0.15. 

superconducting state ui res and Ao yield the kink struc- 
ture along (w, 0) — > (tt, tt) direction that is not present in 
the normal state. 



B. Doping dependence of renormalization 

The different reasons for the kink structures in hole- 
doped cuprates along different directions in the first BZ 
will be also reflected in their doping dependence. So 
far, the results we have shown were for optimal doping 
concentration x = 0.15 that refers to a band filling of 
n = 0.85 31 . Note, the superconducting transition tem- 
perature T c behaves differently in the overdoped (OD) 
and underdoped (UD) regime: 

T c oc A(T -»()), OD 



G 



T c oc n s (T 0), UD 

where n s is the superfluid density that is calculated self- 
consistently from the generalized Eliashberg equations 9 . 

In the antinodal (0, w) — * (ir, it) direction the kink is 
only present below T c due to the feedback of (p(u>). In 
the OD case, 4>(ui) decreases reflecting a mean-field-like 
behavior. Thus, the energy where the kink occurs must 
decrease with overdoping: 



UMnk(x) OC A (x) 



(12) 



This behavior is indeed observed by Dessau and co- 
workers 14 . Note, the above argument remains true also 
in the strongly OD case where no resonance peak in 
Im x(Q,o>) occurs because the feedback effect of (f>(cu) 
should always be present. 

Regarding the kink along the nodal (0, 0) — > (w, n) di- 
rection we note the following. In Fig. 8 we show the cal- 
culated doping dependence of Im \ at the antiferromag- 
netic wave vector Q versus frequency in the normal state. 
One clearly sees the characteristic Ornstein-Zernicke be- 
havior (see Eq. (2)) of Im x, 



Im x(q = Q,w) oc 



U 2 +L> 2 



(13) 



and that u> s f increases with increasing doping from un- 
dcrdoped to overdoped cuprates. Since w s f determines 
the kink position along (0, 0) — ► (ir, ir) direction we ex- 
pect 



Ukink{x) oc uj sf (x). 



(14) 



This is in qualitative agreement with experimental data 11 
(for underdoped regime and optimally doped super- 
conductors). On the other hand, the spectral weight 
of Im x(Q, w) decreases drastically with overdoping. 
Therefore, the coupling of the quasiparticles to spin 
fluctuations is getting much weaker in the OD case. 
These two competing effects seem responsible for the non- 
monotonic and weak doping dependence of the kink po- 
sition in the nodal direction 14 . 

In Fig. 9 we illustrate the kink feature resulting from 
the renormalization (^- ~ v* ~ (1 + A), v ~ v F ) of the 
bare dispersion. We estimate 2 < A < 3. This renormal- 
ization is doping dependent and stronger for underdoped 
hole-doped cuprates. Of course, we expect that the posi- 
tion of the kink as well as the change of the quasiparticle 
velocity (v — > v*) are important fingerprints of the cou- 
pling to spin fluctuations. Note, v* v for frequencies 
uu > LUkink reflects mainly the width of the peak in Im x- 
Important is the slope ratio v* jv for u) < LJkink- 

Another important behavior concerns the asymme- 
try between hole and electron-doped cuprates. Note 
that no kink feature has been reported in the electron- 
doped cuprates 37 . It is believed that the electron-phonon 
coupling is much more pronounced in electron-doped 
cuprates than in hole-doped ones. This is indicated, for 
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FIG. 8. Calculated paramagnon spectrum, i.e. the dynam- 
ical spin susceptibility Im x(Q,w) at a temperature T = 2T C 
for different doping concentrations, x = 0.12 (underdoped), 
x = 0.15 (optimal doping), and x = 0.18, x = 0.22 (over- 
doped). 
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FIG. 9. Illustration of the kink position. We also show 
Re E, since this controls the kink position and v* . The dashed 
curve refers to the bare dispersion. Note that v* — * v for 
u> > LOkink reflects mainly the width of the peak in Im x- 



example, by the behavior of the resistivity p oc T 2 in 
the normal state at optimum doping and by the tran- 
sition between d a .2_ y 2-wave symmetry of the supercon- 
ducting gap towards anisotropic s-wave as it has been 
observed in several experiments 38 . Simply speaking, the 
spin fluctuations in electron-doped cuprates are weaker 
than in the hole-doped ones yielding a smaller T c and a 
smaller superconducting gap 39 . Thus, no kink is present 
in the nodal direction and also no kink occurs in the 
(0,7r) — > (tt,it) direction below T c . This is related to 
the fact that the flat band around (0, ir) lies in electron- 
doped cuprates well below the Fermi level and, therefore, 
it cannot be softened due to <f>(u>). 
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FIG. 10. Calculated feedback of superconductivity on the 
spin susceptibility Im for the electron-doped (a) and 

hole-doped (b) cuprates at optimal doping (x=0.15). The 
solid curves refer to the normal state (T=1.5T C ), while the 
dashed curves denote the renormalized spin susceptibility 
in the superconducting state at T=0.7T C . Due to large 
uj s f = 0.3t and the small superconducting gap, the feedback of 
superconductivity is small in electron-doped cuprates. Con- 
trary, due to a small uj s f — 0.09i in the hole-doped cuprates 
the feedback of superconductivity fulfills a resonance condi- 
tion for Im x yielding a strong renormalization of the spin ex- 
citation spectrum and to a formation of the resonance peak. 
Note that the hopping integral t is different for hole- and elec- 
tron-doped cuprates as discussed in the introduction. 



C. Relation of kink and resonance peak 

In Fig. 10 we show results for the spin susceptibil- 
ity Im xrpa{Q,, in the optimally electron(a)- and 
hole(b)-doped cuprates in the normal and superconduct- 
ing state taking into account different tight-binding en- 
ergy dispersions 19,39 . While in the normal state of hole- 
doped cuprates u) s f is of order of 25 meV, in the electron- 
doped ones its value is much larger (oj s f rts 70 meV) and 
Im x is much less pronounced. Therefore, antiferromag- 
netic spin fluctuations are much weaker in the electron- 
doped cuprates due to weaker nesting of the Fermi surface 
and less density of states due to the flat band well be- 
low the Fermi level. Regarding the superconducting state 
note, that in the hole-doped cuprates a strong renormal- 
ization of the spin fluctuation spectra occurs due to the 



feedback effect of superconductivity and that Ao ~ u> s f 
leading to a resonance peak at w = uo res (see Fig. 10(b)). 
To be more precise, a resonance condition 



1 



= Rc xo(q = Q,w 



U) r 







(15) 



which signals the occurrence of a spin-density-wave col- 
lective mode, must be fulfilled in order to observe a res- 
onance peak 35 . In electron-doped cuprates, the spin ex- 
citations do not obey Eq. (15) and thus only a rear- 
rangement of spectral weight occurs below T c , but no 
resonance peak. Therefore, the kink feature is intimately 
connected with the resonance peak. As we see from Fig. 
10(a) there is only a small feedback of superconductivity 
below T c on Im \ m the electron-doped cuprates due to 
iv s f » A . Thus, we find also no kink feature in the 
superconducting state of electron-doped cuprates in the 
antinodal direction. 



D. Anisotropic scattering rates 

Finally we discuss the anisotropy of the scattering rate 
t~ 1 (uj) of hole-doped cuprates at different points on the 
Fermi surface. In Fig. 11 we show our results for t _1 (w) 
at the antinodal point and the nodal point, respectively, 
for optimal doping (a) and for the overdoped case (b) 
for various temperatures. In Fig. 11(a) one clearly sees 
that the scattering rate is very anisotropic on the Fermi 
surface reflecting the anisotropy of the coupling of ele- 
mentary excitations to spin fluctuations. In particular, 
t _1 (w) in the normal state is almost three times larger 
at the antinodal point than at the nodal point. This 
agrees with recent ARPES experiments 40 . Furthermore, 
we find that Im £ oc lu demonstrating a non-Fermi liq- 
uid behavior in the hole-doped cuprates. Below T c at 
the antinodal point t _1 (cj) reveals a strong feedback of 
superconductivity at energies Lu res + Aq. At the nodal 
point the effect of superconductivity is rather weak. In 
the overdoped cuprates the anisotropy between nodal and 
antinodal points is strongly reduced and for lu — > al- 
most disappeared. Most importantly the system then be- 
haves more Fermi-liquid-like. The latter is seen from Fig. 
1 1 (b) where one observes a crossover from the Im £ oc lu 
to the Im £ oc to 2 behavior. This is also in agreement 



with experimental observation 41 . 



IV. SUMMARY 

In summary, we have analyzed the elementary excita- 
tions in hole- and electron-doped cuprates and the finger- 
prints of spin fluctuations on them. The quasiparticles 
around the antinodal points of the BZ experience the 
strongest scattering on spin fluctuations yielding a non- 
Fermi liquid behavior. In agreement with experimental 
data, we find that coupling of holes to spin fluctuations 
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FIG. 11. Scattering rate t _1 (u) of optimally (a) and over- 
doped (b) hole-doped cuprates versus frequency at the nodal 
and antinodal point of the Brillouin Zone (BZ) calculated at 
various temperatures. The anisotropy results from coupling 
to spin fluctuations and is disappearing in the overdoped case. 
Thus, a crossover from a non-Fermi liquid to a Fermi liquid 
behavior occurs. Note also the feedback effect of supercon- 
ductivity for different parts of the BZ at optimal doping. 



yields a kink feature in the renormalized energy disper- 
sion. 

Possible phonon contribution to the kink fea- 
ture One of the interpretation of the kink structure in 
hole-doped cuprates has been the electron-phonon inter- 
action suggested by Lanzara et al. 13 Indeed, it is clear 
that phonons would also cause a kink structure in the 
energy dispersion if one assumes that Eliashbcrg func- 
tion a 2 F(q, uj) has the same features as x(q, w), namely 
peaked at the wave vector Q and at the Debye frequency 
Wd, i.e. lo = lojj « uj s f. By analyzing Fig. 2 it is 
clear that both spin fluctuations and electron-phonon 
coupling can cause a kink structure. However, in gen- 
eral, one would expect that its position and doping de- 
pendence might be different in both cases. For example, 
only in the case of dominant spin fluctuation coupling 
the kink structure can be related to INS experiments, 
i.e. Im x(Q,w), and, furthermore, the kink position is 
given by Ukink ~ -Ek-Q + uj s f(x). As discussed earlier, 
the kink feature along the antinodal (0, ir) — > (ir, tt) direc- 



tion results from the structure in 4>(u). Thus, additional 
structure in 4>(tv) due to the electron-phonon interaction 
(EPI) may also contribute. Therefore, the question re- 
mains: How to distinguish between spin fluctuations and 
phonons as a reason for the kink formation? To answer 
this question one has to understand how consistent are 
both scenarios with available experimental data. For ex- 
ample, as was shown Zeyher and Greco 36 , the value of 
electron-phonon coupling extracted from the kink anal- 
ysis yields the value of the electron-phonon coupling A 
and T c that are too low to account for high-T c supercon- 
ductivity in the cuprates. Furthermore, assuming that 
the kink structure arises only from the EPI, it is difficult 
to understand the dj.2_ y 2-symmetry of the superconduct- 
ing order parameter and related observed anisotropy of 
the kink structure (see Appendix B). Note, only the 
spin fluctuation scenario yields T c « 70K 30 ' 9 , a d x i_ y i- 
wave order parameter, and a kink feature in qualita- 
tive agreement with experiment. Also the doping de- 
pendence of the kink is difficult to explain within the 
phonon scenario. In contrast to Eq. (10) one would 
expect LOkink ~ -Ek-Q + lod{x) in the case of electron- 
phonon coupling. 

Kink structure in the triplet superconductor 
Sr2RuC>4? Finally we want to emphasize that the 
formation of the kink feature due to spin fluctuations 
should not be restricted to cuprates. For example, the 
quasi- two-dimensional triplet superconductor S^RuO^ 
(isostructural to La2Cu04) 42 reveals pronounced incom- 
mensurate antiferromagnetic spin fluctuations at the 
wave vector = (27r/3, 27r/3) and frequency u) s f « 
6meV that originates from the nesting properties of the 
quasi-one-dimensional a and /3-bands 43-45 (see Fig. 12 
for an illustration) . On general grounds one would expect 
a kink structure in the renormalized energy dispersion of 
the quasiparticles. Although the correlation effects are 
weaker in Sr 2 Ru04 (U is smaller), and Qi is an incom- 
mensurate wave vector, similar conditions as in cuprates 
are present. Note, the kink feature should occurs at 
smaller energies than in cuprates due to a lower value 
of oj s f in the ruthenates. Further experimental studies 
should test our suggestion. 
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FIG. 12. Illustration ol a possible kink structure in 
Sr2Ru04. The Fermi surface of Sr2Ru04 consists of three 
bands. The nesting properties of the /3-band yield to a for- 
mation of two-dimensional incommensurate spin fluctuations 
at Qi = (27r/3, 2-7t/3) and ui s f « 6 meV. Therefore the quasi- 
particles at the /3-band should be strongly renormalized due 
to coupling to spin fluctuations. 



APPENDIX A: GENERALIZED ELIASHBERG 
EQUATIONS 

The interdependence of elementary excitations with 
spin excitations leads to strong self-energy effects. The 
corresponding Dyson equation yields the dressed 2x2 ma- 
trix Green's function G in terms of the bare Green's func- 
tion Go and the self-energy S: 

= iuj n Z(k)T O -[e(k)+t(k)}T 3 -0(k)T 1 , (Al) 

where k = (k, ioj n ). In the FLEX approximation for the 
Hubbard Hamiltonian the sclf-cncrgy £ is determined by 
the following generalized Eliashberg equations: 

E(fc) = 

]T [P„(k - k')T G{k')T + P c (k - k')T 3 G{k')T 3 ] 

k' 

= J2 v Mk-k')G(k'). (A2) 
k' 

In order to provide a better understanding of our nu- 
merical procedure we show the corresponding Fcynman 
diagrams for V e g in Fig. 13. Within RPA the spin and 
charge fluctuation interaction are given by 



P s = (2 7 r)- 1 C/ 2 Im (Sxs-Xso) 
with Xs = XsoO- -Uxso)^ 1 and 

P c = (2ir)- 1 U 2 lm (3 X c-Xco) 



(A3) 



(A4) 



with Xc = Xco(l + Uxco)^ 1 - Therefore, the kernel I and 
the spectral functions of the one-particle Green's function 
in Eq. (5), A v , read 



I{u>, fl, u>') 



/(-a/) + 6(0) f(aj') + b(n) 



uj + iS — CI — lj' lu + iS — Q 



A„(k, uj) = — 7T x Im [a„(k, w)/.D(k, w)] , 



and 



D^^Zf-^ + tf 



a a = luZ, a 3 = e£ + £, ai 



(A5) 



(A6) 



(A7) 



(A8) 



In Eq. (A5) / and b are the Fermi and Bose distribution 
function, respectively. Finally, the bare susceptibility is 
calculated from 

7T f°° 

Im XsO,cO = -T7 / duj' [f(u) - f(u> + u)] 

J —oo 

x ^[7V(k + + w)JV(k, w') 

k 

±Ai(k + q,o/+w)i4i(k,o/)] , (A9) 

where we assume that the same itinerant carriers are re- 
sponsible for the elcmenatry excitations and, at the same 
time, generate the spin excitations. In Eq. (A9) we 
use N(k,w) = A)(k, u) + A 3 (k,uj), and the real parts 
are calculated with the help of the Kramers-Kronig rela- 
tion. The subtracted terms in P s and P c remove a double 
counting that occurs in second order. Note that V e g in 
Eq. (5) is dominated by the exchange of spin fluctua- 
tions due to the fact that the system is in the vicinity 
of an antiferromagnetic phase transition, but the above 
equations still remain valid in the case where Xc becomes 
more important. 

Our numerical calculations are performed on a square 
lattice with 256 x 256 points in the first Brillouin Zone 
and with 200 points on the real w-axis up to I6t on a log- 
arithmic scale. Within our self-consistent procedure the 
full momentum and frequency dependence of quantities 
is kept. 



APPENDIX B: PHONONS AND D x2 Y 2-WAVE 
SUPERCONDUCTIVITY 

In this Appendix we analyse how the magnetic mode 
which is mainly peaked at q = Q = (n, n) leads to 
a d a .2_j / 2-wave order parameter that is maximal around 
(n, 0) and, in particular, to which extend phonons con- 
tribute to this result. In general, the generalized Eliash- 
berg equations read after the inclusion of attractive 
phonons (branch i) via their spectral function a 2 F,(q, ft): 
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FIG. 13. Particle-particle channel ol the Bethe-Salpeter 
equation for superconductivity due to an effective pairing in- 
teraction V e ff entering in Eq. (5). The solid lines refer to G 
and the dashed lines denote the effective Coulomb interac- 
tion U of Eq. (4). Vertex corrections that would yield to a 
renormalized coupling strength U e g (as indicated in Fig. 1) 
are neglected. The summation of the corresponding bubble 
and ladder diagrams is performed up to infinity. While in 
principle it is possible to treat V e ff{x} an d G(k,a>) on differ- 
ent levels, we assume that both quantities are generated by 
the same itinerant quasiparticles. Note, V e s refers to the ex- 
change of spin and charge fluctuations yielding a d a ,2_j / 2-wave 
instability of the normal state. 

£«(k,u;) = 

N- 1 V / dtt Kff (k - k', SI) - a 2 F,(k - k', 0) 

/+oo 
dcj'I{u),n,u)')A v (k',u)') . (Bl) 
-oo 

For a 2 Fi(q, fl) wc employ a Lorentzian in frequency 
around flo « cud (Debye frequency), and a normalized 
form factor i^(q) peaked at q = q pa j r as indicated in 
Fig. 14. The spin fluctuations that are dominating 
V c ff(q,Lo s f) are peaked at q = Q pa i r - 

It is instructive to write down the weak-copuling limit 
of the fi-component of Eq. (Bl) that reads (T = 0) 

4(1) = _ E N^M A(l) , (B2) 

k' k 

where again — ^/A 2 (k) + ef. is the dispersion of the 
quasiparticles in the superconducting state. Note that 
the contribution to the pairing potential is repulsive for 
spin fluctuations and attractive for phonons, respectively. 
In the case where no phonons would contribute to the 
Cooper-pairing (a 2 i^(q) = 0), Kff(q) bridges parts of 
the Fermi surface where the superconducting order pa- 
rameter has opposite signs. This momentum dependence 
of the pairing interaction is indeed required for solving 
Eq. (Bl) and is typical for unconventional supercon- 
ductivity. Note that for a repulsive and momentum- 
independent pairing potential, Kff(q) = const, no so- 
lution of Eq. (Bl) can be obtained. 

How is the kink related to the pairing mechanism? 
Physically speaking, the interdependence of elementary 




FIG. 14. Illustration of d^.^-wave Cooper-pairing for a 
fixed frequency Q — Qo ~ w s / ~ cod due to spin fluctuations 
peaked at momentum k k' = q = Q pa ir and due to phonons 
peaked at q = q pa i r . The solid line denotes the Fermi surface 
and the dashed line refers to the nodes of the d x 2_ y 2-w&ve 
order parameter. The corresponding sign of the order param- 
eter is also displayed. 

excitations that dominate 14s (q), leads to d x 2_ y 2-wave 
Cooper-pairing as well as to the kink structure as ob- 
served by ARPES experiments. In other words, the 
quasiparticles around the hot spots couple strongly to 
spin fluctuations that leads (a) to a 2 -wave order 

parameter, and (b) the same coupling leads to the kink 
in the nodal direction that occurs close to the Fermi level 
where Q pa j r = (tt,tt) as indicated in Fig. 2. 

It follows also from Eq. (14) that attractive phonons 
with a corresponding spectral function a 2 F(q) peaked at 
q = q P ai r contribute constructively to d x 2 _ y 2 -wave pair- 
ing as long as the main pairing interaction is provided by 
spin fluctuations. However, the kink close to the antin- 
odal points occurs only below T c and is a result of 4>{w) 
that is maximal around (0, 7r). Therefore, the kink struc- 
ture in the antinodal direction is mainly connected to 
spin excitations peaked at Q pa ir = (t, 7t) and not to the 
phonon branch peaked at q pa i r . 

Note, in the case where no spin fluctuations would be 
present, i.e. 14fr(q) = 0, the attractive phonon contri- 
bution will cancel the minus sign on the RHS of Eq. 
(Bl) yielding an order parameter with s-wave symme- 
try. Thus, we safely conclude that both, d a .2_ y 2-wave 
Cooper-pairing and the anisotropy of the kink feature in 
the elementary excitations are hardly to recoincile within 
the same physical picture. 
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